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Abstract 

We consider quantum dynamics of the order parameter in the discrete pairing model (Richardson 
model) in thermodynamic equilibrium. The integrable Richardson Hamiltonian is represented as 
a direct sum of Hamiltonians acting in different Hilbert spaces of single-particle and paired /empty 
states. This allows us to factorize the full thermodynamic partition function into a combination 
of simple terms associated with real spins on singly-occupied states and the partition function of 
the quantum Xy-model for Anderson pseudospins associated with the paired/empty states. Using 
coherent-state path-integral, we calculate the effects of superconducting phase fluctuations exactly. 
The contribution of superconducting amplitude fluctuations to the partition function in the broken- 
symmetry phase is shown to follow from the Bogoliubov-de Gennes equations in imaginary time. 
These equations in turn allow several interesting mappings, e.g., they are shown to be in a one-to- 
one correspondence with the one-dimensional Schrodinger equation in supersymmetric Quantum 
Mechanics. However, the most practically useful approach to calculate functional determinants is 
found to be via an analytical continuation of the quantum order parameter to real time, A(r — t- it), 
such that the problem maps onto that of a driven two-level system. The contribution of a particular 
dynamic order parameter, A(r), to the partition function is shown to correspond to the sum of 
the Berry phase and dynamic phase accumulated by the pseudospin. We also examine a family of 
exact solutions for two-level-system dynamics on a class of elliptic functions and suggest a compact 
expression to estimate the functional determinants on such trajectories. The possibility of having 
quantum soliton solutions co-existing with classical BCS mean-field is discussed. 
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I. INTRODUCTION 



The concept of spontaneous symmetry breaking is one of the cornerstones of modern 
physics: Most phase transitions we know are associated with the appearance of a non-zero 
local order parameter that represents a broken symmetry and leads to a state that has a 
lower symmetry than that of the underlying Hamiltonian. In elementary particle physics, the 
Anderson-Higgs mechanism is the most promising scenario to explain the appearance of finite 
masses for elementary particles, including gauge bosons. The canonical model to explain the 
origin of the broken symmetry phenomenon usually involves a Lagrangian for a boson field, 
A, that has quadratic and quartic terms and that can be symbolically represented as follows: 
L[A] = a |A|^ + /3 |A|^ + c |-DA|^, where D corresponds to a gauge-invariant derivative and a, 
/3, and c are constants. In the context of elementary particle physics, it defines a Mexican- 
hat model for the Higgs boson, which is a minimal renormalizable field theory that produces 
symmetry breaking "by design." In solid state physics, such a Lagrangian is associated with 
the Ginzburg-Landau functional for a fluctuating order parameter near a phase transition 
and in many cases it can actually be derived from a more general microscopic Hamiltonian 
(which is typically an interacting fermion model, such that the order-parameter field is 
associated with a composite, rather than canonical boson). 

Such a microscopic derivation was first accomplished by Gor'kov,^ who starting from 
the BCS Hamiltonian obtained the Ginzburg-Landau functional for a superconductor and 
found explicitly the Ginzburg-Landau coefficients in terms of microscopic parameters (z.e., 
electron mass, electron density, interaction strength, and concentration of impurities). The 
general framework for a derivation of this type now appears in excellent textbooks^ and 
can be briefly summarized as follows: One starts with an interacting electron model that 
has a "desired" phase transition (e.g'., electrons with attraction for superconductivity): The 
partition function of the model can be expressed in terms a path integral of the corresponding 
imaginary-time (Grassmann) action, which includes a quartic term describing interactions. 
This term in path integral can be decoupled via an auxiliary Hubbard-Stratonovich boson 
field, A(x) = A(r, r). Then, the fermionic component of the action becomes Gaussian and 
the fermions can be integrated out to produce an effective action S'eff[A(x)], which can be 
formally expressed as a non-linear functional determinant [see, e.g.^ Eq. (12. 6 p in Sec. [H]. 
The Hubbard-Stratonovich field, A(a;), describes a fiuctuating in space and imaginary time, 
r, order parameter and the appearance of a non-zero expectation value. A, of this field 
below a certain transition temperature, Tc, is associated with a broken symmetry phase. 
In the vicinity of Tc, the relevant trajectories of A(a;) are assumed to be such that its 
imaginary-time dependence is unimportant [that is, A(r, r) is assumed to be independent 
of r], A(r) is in some sense small, and it is also assumed to be weakly fiuctuating in space 
(long- wavelength approximation). Hence, the action can be related to the free energy by 
simply writing J-'[A(r)] = TcS* [A(r)], and expanded in a Taylor series, which yields the 
Ginzburg-Landau theory, with the quadratic coefficient a oc (T — Tc). 
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The derivation of the Ginzburg-Landau theory outhned above is justified only near a 
classical phase transition. Below T^, the assumptions about A being small and r-independent 
break down (if the relevant interaction constant, g, is not small they may break down even 
"earlier"). However, it is exactly the low-temperature phase, including the ground state 
that we associate with a spontaneously broken symmetry. This picture is based on the very 
reasonable assumption that the relevant "trajectories" of the order-parameter field, A(a;), at 
low temperatures are located near the classical saddle-point A = const, which becomes the 
only possible trajectory at T = and therefore represents an exact solution. This assumption 
is equivalent to stating that the effective action. See has one and only one minimum 

which occurs in a single "point" in the space of all allowed functions, |A(r, r)| (modulo the 
overall phase). We reiterate that there is no good reason to expect that the simplified form 
of the Ginzburg-Landau action remains reliable at low temperatures. In fact, if we "insist" 
on the canonical Ginzburg-Landau form and attempt to derive the corresponding coefficients 
in the expansion, we shall find that the coefficient of the quartic term generally diverges as 
T — )■ 0.-1^ Hence, we have to work with the full functional determinant in S^s [A(x)], which 
is a complicated non-linear functional and we know little about its properties apart from 
its behavior in a tiny sub-space of constant functions. To the best of author's knowledge, 
there is no model (associated with breaking of continuous symmetry below Tc > 0), where 
such functional determinants have been explicitly calculated beyond the classical mean-field 
analysis. 

The objectives of this work are to bring up the general problem of non-perturbative 
quantum dynamics in broken-symmetry phases and to construct a general framework to 
calculate functional determinants that appear in the non-linear effective action for quantum 
trajectories of the order parameter in the pairing model. The latter is a seemingly hopeless 
goal, but we show that one can obtain exact results in certain cases and based on those results 
formulate a more general Ansatz that is expected to be useful for a large class of quantum 
trajectories. To address this and other related questions, we employ the Richardson pairing 
model,-"- which is an interacting fermion model that has a paired ground state built-in. In 
fact, it is "almost" the mean-field BCS model in the sense that the corresponding order 
parameter does not have any real-space dependence and so all such fiuctuations^ have been 
eliminated. However, the model still retains quantum dynamics of A(r). The Richardson 
model is integrable and there exists an exact Bethe-Ansatz solution,-!^ which determines 
the exact eigenstates and spectrum of the model in sectors with a fixed number of single- 
particle excitations and Cooper pairs. However, this algebraic Bethe Ansatz solution does 
not appear to be very helpful in calculating the thermodynamic partition function in the 
grand-canonical ensemble and we use an alternative method, which is based on coherent- 
state path-integral representation of Anderson pseudospins,- describing the BCS sector of 
the model. We use a mapping of the equilibrium problem in imaginary time onto that of non- 
equilibrium superconductivity, and take advantage of the exact non-equilibrium solutions, 
obtained recently in a series of amazing papers by Levitov et a/.— and Yuzbashyan et alM^ 
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By analyzing a certain family of exact results, we propose a general closed expression to 
estimate the corresponding functional determinant, which is not always exact but is expected 
to be quantitatively reliable for a large class of elliptic functions and their limits. 

Our paper is structured as follows: In Sec. [Tll we present the canonical Richardson model 
and formulate in more technical details the key questions within the conventional Grassmann 
path integral/Hubbard-Stratonovich approach. The questions involve studying various as- 
pects of fluctuation physics and they are addressed in the rest of the manuscript using a 
variety of techniques: In Sec. Illlt we derive combinatorially an exact expression for the 
thermodynamic partition function of a generalized Richardson model in terms of a "spin 
partition function" associated with single-particle states and an "Anderson pseudospin par- 
tition function" associated with the paired/empty states. The generalized Richardson model 
includes the canonical Richardson model (reduced BCS Hamiltonian) as a particular case, 
and in this limit, the spin part of the partition function becomes trivial, so that the problem 
reduces to the problem of calculating contributions of Anderson pseudospins to the parti- 
tion function. Sec. HV] formulates a coherent-state path integral for Anderson pseudospins 
to calculate the functional determinants of interest. It is shown that by introducing a single 
Hubbard-Stratonovich field one can represent the full thermodynamic partition function as 
a product of terms local in parameter space. The contribution of each such local term to the 
partition function follows from the Bogoliubov-de Gennes equation in imaginary time. In 
Sec. lYl we study phase fluctuations within the path integral formalism and obtain an exact 
expression for the partition function in terms of a sum of phase winding numbers. Sec. IVII 
is the main part of the paper, which addresses the question of (possible) fluctuations of 
the amplitude of the order parameter, assuming that the phase fluctuations are completely 
suppressed. Sec. |Vl] contains several parts: In Sec. IVI At the symmetry properties of the 
imaginary-time Bogoliubov-de Gennes equations are discussed and it is shown that the full 
density matrix solution satisfying the proper initial condition, p(r — 0) = 1, can be con- 
structed from a particular spinor solution satisfying arbitrary initial conditions. In Sec. IVI B[ 
we show that the general problem of solving imaginary-time Bogoliubov-de Gennes equa- 
tions in the presence of a quantum-fluctuating order-parameter fleld is equivalent to that of a 
one-dimensional supersymmetric Schrodinger equation, with "superpotentials" determined 
uniquely by A(r). Therefore, the cases where these two problems are solvable are shown 
to be closely related. Sec. IVI C II derives an exact expression for the full density matrix, 
p(r), corresponding to a non-trivial dynamic order parameter, representing the soliton of 
Ref. analytically-continued to imaginary time. The resulting functional determinant is 
found to be surprisingly simple and is equivalent to that of a Fermi gas. Sec. IVI CI suggests 
that the simpliflcation of the functional determinant observed in Sec. IVI CTl is not accidental 
but has a natural explanation: It is argued that the effective action associated with a given 
quantum-fluctuating, A(r), is given by the sum of the dynamical phase and Berry phase 
accumulated by a two-level-system driven by a time- dependent magnetic fleld determined 
by the analytically continued order parameter, A(r — it). This conjecture is verifled to 
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work well on a large class of functions, where A(r + it) is an elliptic function with two 
primitive periods along the r and zt-axes. A general expression for the corresponding effec- 
tive action is presented in Sec. IVI Dl and the possible implications of the results obtained to 
non-perturbative quantum dynamics of the superconducting order parameter are discussed. 



II. THE RICHARDSON PAIRING MODEL AND KEY QUESTIONS 



Let us consider spin- 1/2 fermions, described by the creation/annihilation operators, c] ^ 
and ci^s, labeled by the spin index s = ±1 and the index / G £, where £ is a set of allowed 
single-particle states. It can be a discrete, possibly finite, set (associated for example with 
localized levels in a mesoscopic superconducting grain^^"— ) or a continuum of momentum 
states in a system with open boundary conditions (such that |/, s) and |/, —s) are a pair of 
time-reversed states). We will refer to the states / as to "sites." We perform some formal 
mathematical manipulations assuming that C is discrete and finite, but it is without a loss 
of generality, as this assumption does not preclude us from taking the proper limit at any 
stage of the calculation. The canonical Richardson Hamiltonian (or equivalently the reduced 
BCS Hamiltonian) describing an s-wave superconductor has the form: 



iec,s=± 



9 
2Vc 



(2.1) 



where Vc is either the number of sites in C, if the set C is discrete, or otherwise if C represents 
a continuum spectrum, Vc is a volume (in this case, the sums are to be replaced with integrals 
over momenta, k = /). In what follows, we will also use the notation g = g/ {'^Vc). 

To formulate the main questions, let us first follow the conventional method of treating 
Hamiltonian (12. ip and represent the partition function as a Grassmann path integral: 



/3 



Zr= I Vci^sV\Dci^s[r]QW 



i,i'ec 



i^i - dr) ci^s - gci,+ci-cv-cv,+ > . (2.2) 



and introduce 



Now, we use the identity, e^'^'^'^'^ = J ^exp -| |A| + (Acc' + Ac'c) 
the Hubbard-Stratonovich field, A(r), to decouple the interaction term in the Grassmann 
action, integrate out the fermions from the resulting quadratic theory, and arrive to the 
following standard effective action expressed in terms of the order parameter field 



Zr= / I?A*[r]I?A[r]e 



JjDet [9^-^;f^-ReAf^ + ImAfJ^ 



(2.3) 



where r are two-by-two Pauli matrices in the Nambu space and the determinant is to be 
evaluated over both the time variable and the Nambu space. To trace over the Nambu space, 
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one can use the identity below 



Det 



A 


B 


C 


D 



det [A ■ D] det [l - D'^ ■ C ■ A'^ ■ B] 



(2.4) 



which is vahd for any matrices/operators A, B, C, and D, provided that A and D are 
invertible. Applying this identity to Eq. (I2.3p . we find 



Zr = Zpg J VA*[t]VA[t] expj-iy" rfr | A(r)| V 5efr [A(r)] 

where the effective action is given by 

^eff [A(r)] = 5^Tr In [1 - G+ ■ A* ■ G; ■ A] 



(2.5) 



(2.6) 



lec 



and ZpG = Yliec l^r " Cf] is the partition function of a non-interacting Fermi gas given 
by Eq. f l3.18p below. Here Gf = {dr ^ are Green functions, whose explicit form in 
r-representation is easy to obtain. 

Calculating formally the first variation of the effective action 5*0^ [A] with respect to A[r] 
leads to the mean-field equation for an extremum Amf[t] of the functional, which generally 
has the complicated operator form: 

4A^p(r)5(r - r') = 5^ [1 - G+ • A^p ■ G; ■ Amf]"' • [G+ ■ A^p • G;] ir,r'), (2.7) 
^ lec 

where the right-hand-side is to be understood as a kernel of the corresponding operator in 
i.e., a kernel, K{t,t'), defines an operator by its action on an arbitrary 



The equation fl2.7p 



r-representation 

/3-periodic function, /(r), as follows K ■ /(r) = K{T,T')f{T')dT' 
can be cast into a more friendly form of an integral equation, but it would still remain too 
complicated for a systematic analysis. We do however know that there exists a solution 
to this equation, which is a constant that in the classical BCS model is given by Abcs ~ 



Woe 



-1/(^^9) 



(here we have to assume that C is momentum space and ^ ■ = V^z/ J d^r, v is 



lec 



the density of states at the Fermi level, and is the usual high-energy cut-off to regularize 
the Cooper logarithm). One can verify explicitly that indeed Abcs is a true minimum [z.e., 
it is not only a minimum on a tiny subset of constant functions, but a true minimum on the 
space of allowed functions, A(r) = A(r + /3)], but there are still a few important questions 
that remain: (i) Does the classical BCS mean-field result represent the only minimum at 
< T < Tc, or there may exist quantum non-perturbative trajectories of A(r), which 
would give contributions energetically comparable to the classical mean- field (or better)? 
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(ii) A related key technical question is whether it is possible to calculate the functional 
determinant, det [l — ■ A* ■ Gj ■ A] , for "trajectories" of the order parameter with non- 
trivial quantum dynamics? (iii) What are the effects of quantum fluctuationsi^ii^ of the 
modulus and/or phase of the order parameter on thermodynamics [e.g., the energy of the 
ground state)? We will address these questions to some extent in the following sections 
using an alternative method, namely the path-integral formalism for Anderson pseudospins. 

III. FACTORIZATION OF THE GENERALIZED RICHARDSON HAMILTO- 
NIAN 

The Richardson Hamiltonian (12.1 p is known to be integrable^'^ and its integrability is due 
to the existence of an infinite number of conservation laws at two levels of the problem: First, 
the Hamiltonian commutes with the z-component of the spin on any site and therefore the 
Hilbert sub-spaces associated with the singly-occupied states and the paired/empty states 
are separated and can be studied independently.— After this factorization, the Hamiltonian 
for paired states reduces to a pseudospin Hamiltonian (expressed in terms of Anderson pseu- 
dospins). As Richardon discovered, the pseudospin Hamiltonian amazingly has an infinite 
number of conservation laws as well and this allowed him to construct an exact Bethe-Ansatz 
solution to the corresponding spin problem in a given sector (with a fixed total pseudospin), 
and in particular, find a set of coupled algebraic equations determining the energy spec- 
trum in the sector. The Richardson equations are exact and therefore include correctly all 
quantum fluctuation effects, but this exactness also makes it difficult to use the solution 
for practical purposes and to interpret its physical meaning, because the solution mixes up 
fluctuations of the order parameter of different types. In addition, the Richardson equations 
are still too complicated to allow a further analytic treatment and most importantly they 
address different pseudospin sectors independently. For this reason, we do not use the results 
of the algebraic Bethe-Ansatz approach to calculate thermodynamic properties of the model, 
but we find it however very useful to perform the first simpler step in the Richardson's so- 
lution, i.e., to factorize the Hilbert space into single-particle and paired/empty states. It 
turns out that this factorization is allowed for a more general Hamiltonian than (12. ip and in 
the interest of generality and future work, we present this procedure for such a more general 
model, which we dub the generalized Richardson model [e.g. Eq. (13. 4p below represents a 
generalized Ising- Richardson model]. 

Let us define the density, spin, and Cooper pair operators on each site as follows: 

Pi = '^^c] gCi^s is the density operator, (3.1) 

s=± 

^i = l^cls^s,s'Ci,^' (3.2) 

^ s=± 
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is the spin, with •S'f = | ^ sc|s^',s being its z-component, and 



s=± 



P/ = - s cJsc|-s is Cooper pair operator. 



(3.3) 



s=± 



Clearly (p/)^ = q,_q,+. 

Let us now use these operators to express the following generalized Ising-Richardson 
model: 



(3.4) 



h,l2&C 



where describes single-particle energy eigenvalues/spectrum, Bi is an applied magnetic 
field in the z-direction, g is an interaction in the BCS channel, and J is an Ising-type spin 
interaction. We reiterate that the special case of Hamiltonian (13. 4 p with J = B = and 
gi^ i^ = g = const yields the canonical s-wave Richardson pairing model (12. ip that we actually 
will study in the rest of the paper. 

However, the more general Hamiltonian (13. 4p has the same "local" in C conservation 
laws, since it commutes with the z-component of the spin, S[, on any site: 



0, VZ e C. 



(3.5) 



This allows us to define the following projectors for an arbitrary subset of C 



A = n (2^f 



leCi 



and 



ieC2 



(3.6) 



(3.7) 



where £i_2 C C. Note also that ['^1,2] = A, 2 ['^1,2]- By convention we shall denote 
the projectors on a single site {i.e., if the corresponding subset consists of a single element, 

£12 = {/}) as follows: A(0 = (25f) and A (A;) = i - (^25f) . Obviously for those 
single-site projectors we have: 

A(/) + A(0 = i- (3.8) 



This resolution of unity allows us to represent the Hamiltonian (13.41) as a sum of Hamiltonians 
acting in different "sectors" of the Hilbert space as follows: 



lec 



CiuC2=c 



(3.9) 



Each term in the above sum represents two Hamiltonians acting on single-particle states in 
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£i and paired/empty states in £2- The corresponding spin and pairing Hamiltonians are 




K^n[Ci]=V,[C,]\Y,UiPi-BiSf)+ Yl JiuiA^U (3.10) 



and 



T^BCS' [A] = ^2 [A] I 5^ - Yl ~9h,i.PtA\- (3-11) 

Now, one can follow Anderson and check that the operators P/, Pi, and (p/ — 1), when 
constrained by the projector on empty/paired states, form a closed su(2) algebra on each 
site (here and below, we use the symbol, su(2), for the Lie algebra and SU{2) for the Lie 
group) or in other words, the operators are Anderson pseudospins. One can therefore drop 
the projectors and replace the operators with Pauli matrices (since, = 0, we have to use 
the two-dimensional representation) P/ = f^^. Pi = , and pi = + 1. Similarly, one can 
remove the projector in Eq. (13.101) and simply replace pi with one, since each site in Ci is 
guaranteed to be singly-occupied by construction. This leads to the following decomposition 
of the Hamiltonian (13. 4p 

nGK= Yl {^spin[/:i]+HBCS'[>^:2]}, (3.12) 



CiUC2=C 



where T^spin [-^i] and T^bcs' [-^2] are spin-1/2 Hamiltonians acting in different Hilbert spaces. 
These spin Hamiltonians are of Ising and XF-type correspondingly: 

HspiA^i] = Y (^^ - l^^^i) + \ E -JhM^h^L (3.13) 

and 

T^BCS' [/:2] = 5^ 6 (rf + 1) - Y ~9h,hK%- (3-14) 

Since the Hamiltonian (13. 4p does not have operators that connect different partitions of 
C, the total partition function is given by a combination of the products of the partition 
functions corresponding to the Ising and XF-models on different sets 

Zgr= Y Z,pin[jCi]x ZbcS'[C2]- (3.15) 
CiuC2=c 

Note that factorization of the Hilbert space into single-particle and pair/empty states, which 
led us to Eq. (I3.15p . does not require that Sf is locally conserved, but requires only that 
the spin and pseudospin sectors can be uncoupled via projectors (13. 6 p and (13. 7p . which is a 
much weaker requirement. This implies that this construction may be applied to even more 
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general Hamiltonians of type (13. 4p . which include quantum interaction terms for real spin. 
This avenue will be explored elsewhere,— but here we instead focus on the much simpler 
canonical Richardson pairing Hamiltonian (12. ip . where there are no interactions for real 
spins {Ji^.i2 = 0)5 '^'^^ there are magnetic fields {Bi = 0), and hence the partition function 
associated with the single particle states is simply Zgpin = Hie/:! (2e~'^^'), so that the 
full partition function of the pairing model is simplified to 

Zr= Yl U [^e-^^''] Zbcs'[C2], (3.16) 

where Z-qcs' ['^2] is the partition function of the XF- model given by T^bcs' [^2] in Eq- (I3.14p 
on a subset C,2 d C- and where gi^^i^ =5'/ (SV^). We will use this decomposition (I3.16P in 
the remainder of the paper. 

To run a simple sanity check on the result obtained, we consider the non-interacting case 
with (7 = 0, z.e., the Fermi gas. Eq. (I3.14p therefore is the Hamiltonian of non-interacting 
pseudospins in magnetic fields, b/ = (0,0,^^), and the partition function is given by 

^FG= E n [^e-'^^'^] X n K^^'^2cosh(/36J]. (3.17) 

Since the partition function involves products of "local" in L terms and all possible decom- 
positions are to be considered, we can equivalently rewrite Eq. (I3.17p as follows 

Zfg = n [2e-''«' + cosh (/^^O] = H + ' ' (^.18) 

lac i&c 

which is indeed the partition function of a non-interacting Fermi gas of spin-1/2 particles. 



IV. PATH INTEGRAL FOR ANDERSON PSEUDOSPINS 

In Sec. Illlt we showed that the full partition function of the Richardson model is given 

by 

Zr = exp[-/3 5^e/] ( n 2) X ZBCs[C2i (4.1) 

where Zbcs is the partition function of the XF-Hamiltonian with infinite-range interactions 
[here, we subtract a constant from the Hamiltonian "Hbcs' given by Eq. (I3.14p and set 
9h,h = 9]- 

T^Bcs [A] = E ^ir[ - ^ ( E X ( E ) • ^^-^^ 
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To calculate the partition function, we employ the coherent-state spin path integral formal- 
ism and write it in the form: 



^BCS [-^2 



[^2 









f ^ 




/ 




exp < 












I 


.ieC2 



(4.3) 

where = (sin 6i cos (pi, sin 6i sin (pi, cos 6i) is a vector constrained to move on a unit sphere, 
rji = cosOi, nf = {nf ±irvl) /2, and dVtiiji) = d(f)i{Ti)drji{Ti) for any imaginary time, Tj G 
[0,/3 = l/T]. 

We now perform the Hubbard- Stratonovich decoupling for the interaction term in the 
spin path integral, which allows us to write the full partition function in the form: 



Zr = exp 



Mr) 
Tig 



e \[{2 + zi[Mr)]}, 

lac 



(4.4) 



where the zi is a "local" path-integral, which depends on a realization of the "global" 
Hubbard-Stratonovich field 



Zi[A{t)] = / Vn{r)exp 



dr 



--(f) + iiri + A(r)n" + A*(r)n^ 



(4.5) 



Note that in Eq. (14. 4p the explicit factorization of the terms into single-particle and 
paired/empty states is no longer necessary due to "locality" of the "dynamic partition func- 
tion," zi, after the Hubbard-Stratonovich decomposition. The contribution of the single- 
particle terms is simply given by the factor of two in Eq. (14. 4p . 

To treat the path integral (14.51) . we note that it can be "generated" as a solution to the 
following differential equation for a "density matrix," pf 



dpi 
dr 



-hi{T)pi 



A*(r) -6 



\pi, with p(0) 



(4.6) 



The trace of the two-by-two "density matrix" evaluated at r = /3 gives the desired partition 
function 

zi = 1ipi{P). (4.7) 

This relation can be proven by writing a formal solution to Eq. (14. 6 p as a r-ordered expo- 
nential and then expressing it as a path integral to reproduce exactly (14. 5p . 

To verify that the formulas obtained so far are consistent with what is known, let us 
consider the case of the classical mean-field, where the order parameter is taken to be a 
constant ABcsMF('i") = Abcs = const. In this case the solution to Eq. (14. 6 p is given by 
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pf^ = exp (^—hifij . Since, hi = + Re A f"^ — Im A f^, and we can write 

pf"^ = i cosh {EiP) - {m ■ t) sinh (EiP) , (4.8) 



where hi = Ei (n^ ■ r) with In^l = 1, so that Ei = y + {A]"^ is the famihar quasiparti- 
cle spectrum in BCS theory, which in the pseudospin language translates into an effective 
magnetic field experienced by a pseudospin. Calculating the trace, we recover the partition 
function of a spin-1/2 in a magnetic field of mag nitude |b;| = Er. zf^ = 2cosh{Ei/3). Now 
returning to Eq. ( 14. 4p and noticing that {2 + 2 cosh {Ei/3)} = {2 cosh } , we can write 
the classical mean-field contribution to the partition function as follows: 



r A fcPA { /3|A|^ r , fEil3\ 

■^BC?,M¥ = (iWy-P 2^Ki\ j — ^exp<^ h2 2_^ln 2 cosh I— -I 



— 12 

/ ^;, H \ I 

(4.9) 



where we recall that g = g/{2Vc). Varying the action with respect to A, we indeed recover 
the familiar BCS self-consistency equation 



1 = 2Vc^E[^ tanh 



Ei(3' 



(4.10) 



We note that eventhough the classic BCS equation follows from the Richardson Hamiltonian, 
this zero-dimensional model does not have a true (classical) phase transition. In particular, 
if we calculate the Riemann integral over A that appears within the classical mean-field 
approximation in Eq. (14. 9p . the resulting function 2'bcsmf(7') will be continuous in the 
vicinity of a nominal Tc {e.g., one can expand the free energy into a Taylor series and 
obtain a zero- dimensional Landau theory, which leads to a continuous partition function 



expressed in terms of the error function; see Ref. 16| for details). If the underlying physical 
model is higher-dimensional, then a phase transition is anticipated,— and we can interpret 
the temperature at which a derivative of the partition function over T has the sharpest 
slope as a temperature where the phase transition occurs. However, it is not only the 
partition function itself that is of primary interest, but also the trajectories that provide 
main contributions to it. In the weak-coupling limit, the transition point can in turn be 
identified (in the leading approximation with respect to g) with the point where a non-trivial 
solution to the self-consistency equation ( I4.10p first appears, but in strong coupling this is 
not necessarily so. We note that one can use the simple BCS result (14. 9 p for estimates of 
Tc{g) by examining the partition function as explained above {i.e., looking for a temperature 
where the slope of its second derivative is the sharpest). However, of course this procedure 
is not quantitatively reliable as it neglects superconducting fluctuations in real space (which 
are classical fluctuations for the purpose of determining T^.), which have been excluded from 
Richardson model from the outset. In what follows, we will not address the very interesting 
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question of determining in strong coupling, but instead will focus on the effects of quantum 
dynamics of the order parameter. 



V. PHASE FLUCTUATIONS 

Let us now express the order parameter in Eq. ( 14. 4p explicitly as a product of a time- 
dependent amplitude part and a phase factor, A(r) = Ao(r)e*'^(^\ To proceed further, we 
note that the first term (the factor of two) in the product in Eq. (14. 4 p originates from single- 
particle states, which are free (real) spins and as such this factor of two is nothing but a 
partition function of a free spin-1/2. In the path-integral language, it can be "generated" 
by the action, which contains just a Wess-Zumino term and no Hamiltonian. I.e., we can 
use the following "representation of the factor of two:" 



Vn{T)e 



(5.1) 



Such Wess-Zumino terms appear in the factors zi in Eqs. ( 14.41) and ( 14. 5 p as well and we get 



Zr = e 



A(r) 
■ng 



.lac 



exp 



dr 



lac 



n 1+ 



lac 



exp 



dr 



iiTji - 2Ao cos (0/ - 7) sjl-ril 



(5.2) 



We note here that the first term in the product, which is equal to one within our conventional 
Richardson model (12. ip . will have a more complicated form in the generalized Richardson 
Hamiltonian (13.41) . where it should be related to the partition function of an Ising model for 
real spins on singly-occupied sites. 

We now perform the following change of variables ("gauge transformation") — )■ (/)i(r) -|- 
7(r). The dependence of the action on the overall phase of the order parameter disappears 
from the last term in the product in Eq. ( 15. 2p and appears only in the Wess-Zumino term. 
The corresponding 7-dependent part of the action therefore reads 



P 

I I drJ^iiT) [l + vi{r)]. 

n lec 



(5.3) 



We can now evaluate the path integral over 7(t), following Ref. 20| and keeping in mind 
the periodic boundary conditions for A(r) and n(r), so that 7(/3) — 7(0) = 27iq, with q E Z. 
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Therefore, we obtain 



f V^{r) e-^- = J] exp i 2z7rg J] [1 + r/,(0)] /2 i 6 



(5.4) 



Hence, the phase fluctuations of the order parameter constrain the sum fliij) to be equal 
to a constant at all times. The resulting sum over q can be rewritten as an inverse discrete 
Fourier transform 



Therefore, the result of path integration in (15. 4p is 



1^76 



(5.5) 



The partition function for Anderson pseudospins (15. 2p reads 



e > 

M=0 



V 



L/e£ 



lec 



-5AQ-5'wZ-'S'efl[{n;}] g-j 



where we limited the sum over Af to positive values only because (1 + 77^) > (if the set C 
is finite we can restrict the sum to A/" < Vc) and the path integral over the order parameter 
field includes only the dynamics of the modulus. In Eq. (15. 6p . S'aq is the "bare action" for 
the order parameter field, S'wz is the sum of all Wess-Zumino terms for the pseudospins, 
and the interacting part of the effective action reads: 



-5eff[{n,}] 



JJ h + exp 



lec 




^ir]i - 2Ao cos(f)i\ 1 - 7]l 



(5.7) 



We see that the effect of phase fluctuations of A(r) is to separate the partition function 
into "sectors," where the total projection of the z-component of Anderson pseudospins is a 
constant integer at all times. This analogy can be made more explicit, if we imagine the 
associated real-time pseudospin dynamics, governed by the Bloch equation, = b; x M;, 
and where the effective magnetic field is determined by b; = (Re Ao(zt), Im Ao(— it), ^i) 
[where Ao(— it) is the modulus of the order parameter properly analytically-continued to real 
times, T ^ it]. The 5-functions in Eq. (15.61) demand that the real-time dynamics of individual 
pseudospins must be correlated in such a way that they pin the "total pseudospin moment," 
^jMf(t), to a constant. Note that these constraints imposed by the phase fluctuations 
are in addition to the constraint that may be imposed by any mean-field treatment of the 
remaining path integral over the amplitude Aq. From this, one can see that our ability or 
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lack thereof to satisfy a certain mean- field (in a mesoscopic integrable systen>i^ii^) , e.g. a 
constant amplitude such as in classic BCS mean-field, is determined by the initial conditions 
for the pseudospins. 

The physical meaning of all these results can be clarified if we first consider the subset 
of paired/empty states and recall that the density operator on a site, /, (13 .ip of the original 
model is given by = 1 + ff for the paired/empty states / G £2 C £, so that the Anderson 
"spin-up" corresponds to the existence of a Cooper pair and a "spin-down" to an empty 
site. Therefore the operator corresponding to the total number of Cooper pairs is given by 

■^CP = ^ 5^ [1 + rf ] (5.8) 

and the time-dependent field in the path-integral formalism corresponding to this operator 
is given by 

^cAr) = ^ E [1 + ^°^^'(^)] = ^ E [1 + ^'(^)] • (5-9) 

From Eqs. (15. 3 P and (15.91) . we see that the action that includes the phase of the order 
parameter can be written as follows 

S,[C2] = -l ! dTj{T)Afcp{T). (5.10) 



I.e., we recover the fact that in the absence of gapless excitations, the phase of a Bose field 
operator and the number of bosons (Cooper pairs in our case) are canonically conjugate 
operators, satisfying therefore the Heisenberg uncertainty principle. However, the effective 
action (15. 7p may contain contributions from single-particle states as well [they are associated 
with the factor of one in the sum in Eq. (15. 7p ] and if we allow such states, i.e., if / G £1 7^ 
0, then the canonical conjugate to the phase, 7, the way it is defined above, will also 
have a contribution from the real spins on singly-occupied sites. The meaning of the field 
[1 + rjii^r)] /2 is different for those singly-occupied states and relates to the ^-component 
of the actual magnetic moment of a site. This suggests an interesting relation for the full 
phase action (15. 3p . which now includes contributions from both paired/empty states and 
single-particle states: 



p 

S^[C] = —i J dT'j{r) 



-M,t(r) + 5, 



2 



tot 



(5.11) 



where A/tot (t) and S^^^{t) are fields corresponding to the total number of particles and the 
total magnetic moment of the system. If single-particle states are completely gapped out 
as it is usually assumed, then all particles are bound in Cooper pairs, the total magnetic 
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moment is identically zero and we recover the familiar conclusion summarized by Eq. (15.101) . 
But in general, the Hamiltonian version of Eq. (15. lip will be a Heisenberg uncertainty 
relation/commutator, which involves both the superconducting part (Anderson pseudospins) 
and a magnetic part (real spins): 7, \N'tot + S^^^ = il. We reiterate here that while our 
model is "biased" towards a superconducting state and has no magnetic interactions for 
real spins, a more general Richardson Hamiltonian [see, e.g., Eq. (13. 4p ] may have non-trivial 
magnetic interactions [see, e.g., Eq. (I3.13p ]. which in principle may lead to a magnetic phase 



transition that would compete with superconductivity, c.f., Refs. [2l|, [22|, [23|, and \24\. 

Both A/tot and S^^^ are certainly good quantum numbers and are separately conserved. 
Hence, the phase, 7 fluctuates strongly (via the Heisenberg uncertainty principle) and since 
we treated these fluctuations exactly in Eq. (15.60 . the 5-function constraints there effectively 
enforce these underlying global conservation laws. An important question is whether we 
actually need to enforce them to describe a realistic superconductor. The classic description 
of an s-wave superconducting ground state requires no gapless excitations (£1 = 0) and 
hence the phase 7 is identified with the phase of a Cooper-pair superfiuid with broken gauge 
symmetry (that is, 7 does not fluctuate). Per the same Heisenberg uncertainty principle, we 
must require then that either A/cp or S^^^ or both fluctuate strongly (in a closed system, it 
must be both, because the only way by which A/cp can change is by breaking Cooper pairs 
into single-particle excitations). 

Another more technical way to argue in favor of the same conclusion is to consider a 
Richardson model or a more general (non-integrable) physical Hamiltonian from which it 
descends, weakly coupled to a bath and/or to a noisy magnetic field. Then, we are allowed 
to break weakly some constraints associated with the global conservation laws. This can be 
done by "softening" the (5-functions in Eq. (15. 6p : E.g., we can represent each 5- function as a 
narrow Gaussian and then allow a finite width to the Gaussian, which would be equivalent 
to introducing a charging-energy-like term to the action 63^ oc J 'j'^dr that penalizes phase 
fluctuations. Both these arguments suggest that to describe a realistic superconductor in 
the actual broken-symmetry phase, we have to suppress phase fluctuations, which can be 
accomplished by dropping the S'^-term and the resulting constraints in the partition function 
(15. 6p . This however brings up the question of whether the low-temperature state with broken 
gauge symmetry will allow fluctuations of the amplitude of the order parameter and if yes, 
whether they are purely mesoscopic or may involve more serious non-perturbative solutions. 

VI. AMPLITUDE FLUCTUATIONS 

A. Bogoliubov-de Gennes Equations in Imaginary Time 

We now consider the amplitude fluctuations assuming that the phase fluctuations are 
suppressed. As it was shown in Sec. IIVI the partition function, originating from a 
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non-trivial fluctuating order parameter fleld, is given by the trace of the density matrix 
2i[Ao(r)] = Tr /)(/?) , which now is the solution to the following Bogoliubov-deGennes equa- 
tion in imaginary time with a real, but generally time-dependent Aq{t), 



dpi 
dr 



Ao(r) 
Ao(r) -ii 



) 



pi, with p(0) = 1. 



(6.1) 



What is required at this stage is to flnd a general expression for zi [Ao(r)] as a functional of the 
order parameter and to perform a variational analysis on the resulting effective action. This 
is equivalent to calculating the functional determinant in Eq. (12. 5p . which appears within a 
more conventional treatment. This is a difficult problem, which is intimately related to the 
problem of dynamics of a two-level system in a time-dependent magnetic fleld (generalized 
Landau-Zener problem), described via non-linear differential equations that have known 
analytic solutions only in a few special cases. While to determine the exact dynamics of 
pseudospins under an arbitrary perturbation, Ao(r), may not be possible, one can still get 
further insights by taking advantage of the recent progress in understanding non-equilibrium 
BCS super conduct ivityi^ and the problem of dissipation due to externally driven two-level 
systems,— where exact solutions can be obtained for a wide class of external perturbations 
associated with elliptic functions. Below, we explore solutions to Eq. (16. ip in some special 
cases and generalize the results to express the functional determinant that arises within this 
class of dependencies in a compact form. 

However, let us start with a general analysis of the imaginary-time Bogoliubov-de Gennes 
equations (16. ip . Let us assume flrst that Ao(r) = Ao(— r), i.e. that it is an even function, 
which may occur "naturally" or via a periodic continuation from the physical imaginary-time 
interval, [0, f3] [all conclusions below can be generalized easily to the case where Ao(r) = 
Ao(2ro — r)]. Let us also consider a Nambu spinor x(r) = (^|^|) and look for a solution to 
the following equations 



without specifying initial conditions. We also require that Ao(0) = Ao(/3), since it is a fleld 
that arises from a path integral in imaginary time. The corresponding function may have a 
"natural" period commensurate with /3 or a single "accidental" period and in the latter case 
we shall periodically continue the function, Ao(r) deflned on r G [0,/3], such that it satisfles 



immediately implies that X2(t) = {—ir^) Xi(— t) (if^) is also a linearly-independent solution. 



drU = iu + Ao(r)f ; 
drV = -^v + Ao(r)-u, 



(6.2) 
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Hence a general solution to Eq. (16. 2p has the form 



where C± are arbitrary constants determined by the initial conditions. Note that to deter- 
mine the "density matrix," pi, in Eq. (16. ip . we need to find two solutions that satisfy the 
initial conditions, iI>i{t — > 0) = (J) and il>2iT — )• 0) = (°). Let u{0) = uq and v{0) = vq 
be the initial conditions of a solution, that we assume known. Then, from Eq. (16. 3p . 

we determine the solution that satisfies the first required initial condition {i.e., spin-up at 
r = 0) as follows 

Ul + Vi \UoviT) - Voui-T)J 

Per the same argument as above, the time- reversed to this solution, iI'2{t) = 
{—ify) iIji{—t) {ify) satisfies the other initial condition {i.e., spin-down at r = 0). Therefore, 
we conclude that if we know any solution to Eq. (16. 2p . we can construct the 2x2 "density 
matrix" as follows: 

p(r) = {^i(r); (-^f^) ^i(-r) (^f^)} , (6.5) 

where the solution, V'i(r), and its time- reversed form the columns of p{t). This results in 
the following "partition function" (functional determinant) of interest: 

These results can be readily generalized to the case, where the order parameter is an even 
function with respect to an arbitrary tq G [0,/3], i.e., if Ao(r) = Ao(2ro — r). The time 
reversal operation that generates another linearly independent solution can be written as 
follows \&2(t, To) = (— ?^^) ^i(— T, — To) (ir^), where \l'i(r, ro) is a solution satisfying the 
initial condition \E'i(0,ro) = (q) and, which itself can be constructed out of an arbitrary 
solution as follows ^{r) = C7+(ro) I-J) + C^-(^o)(7(^-"^)- 



B. Bogoliubov-de Gennes Equations and Supersymmetric Quantum Mechanics 

We see that if we know any particular solution to Eq. (16. 2p with arbitrary initial con- 
ditions, the problem of calculating the functional determinant is solved. However, it is of 
course the main challenge to find a particular solution. To shed light on the complexity 
of this general problem and to obtain a further interesting insight, we now take a detour 
to point out a direct connection between the solvability of Bogoliubov-de Gennes equations 
(16. 2 p and supersymmetric Quantum Mechanics. 
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Let us introduce the following functions 



p{t) = u{t)v{t), i?+(r) = and RJt) - ""^^^ 



u{t) 



v{r) 



From Eqs. f l6.2p . we find 



< d^R+{T) = -2^R+{t) + Ao{t) [1 - Rli, 



(6.7) 



(6.8) 



, drR-ir) = 2^R_{t) + Ao(r) [l - Rl{r)] . 
The function, p(r), is expressed in terms of the other two related functions R+{t)R_{t) = 1: 



p(r) =poexp <^ / dsAo{s) [R+{s) + R_{s)] 



(6.9) 



where po is a constant of integration that can be set to one, po = 1, since we are looking for 
an arbitrary solution. We see that Eqs. (16. 8p for R±{t) are represented by a rather general 
Riccati equation, which has been studied for some 300 years and which has known analytic 
solution only in a limited number of cases. But let us however proceed further and simplify 
the form of these equations by introducing the following new variables: 



x{T} = I Ao{s)ds, W{t) 

'0 



An(r) 



and r±(r) = i?±(r) ± W^(r). 



(6.10) 



We now assume also that Ao(r) does not change sign (which in fact is a requirement if phase 
fluctuations have been eliminated, since a change-of-sign in the order parameter should be 
incorporated into its phase dynamics). In this case, we can unambiguously determine t{x) 
and treat all functions involved as functions of x. We get 



r'^{x) + r|(x) = 1 + W\x) ± W'{x) = 1 + K 



(6.11) 



These are Riccati-type equations as well, but they now have a form that is reminiscent 
to equations appearing in the context of supersymmetric Quantum Mechanics. To see the 
connection, we recall that a generic Riccati equation can always be reduced to the following 
form 



(6.12) 



cr(x) (j){x) ' 

such that a particular solution to Eq. (16.12^ is written explicitly as fo{x) = a / (j) + 



const + J -^dx 



and therefore the question of finding an analytic solution to a 
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generic Riccati equation f l6.12p reduces to that of finding explicitly the functions and 
(7{x). In our case (16.111) . we see that <J±{x) = (p'^{x), while the equations for (f)±{x) have the 
form 



2 



dx' 



+ VM 



^±{x) = -M^). a;e[0,L], (6.13) 



where L = J Ao(s)(is is the period of the potentials, V±{x) = W^{x) ± W'{x), and W{x) is 



defined in Eq. (16.10p . We see that the operators H± in the right-hand-side of Eq. (16.13^ are 
Schrodinger operators associated with two superpotentials V±{x), which have the canonical 
form of those in supersymmetric Quantum Mechanics^ and which in our case are determined 
by the underlying dynamics of the order parameter! Furthermore, since the Ao(t) has 
been periodically continued, Eqs. (I6.13P are actually Schrodinger equations in a periodic 
superpotential. 

Even though, for our purposes what is really needed is the "wave-function" associated 
with just one (negative-energy) state, E = —1, we can easily examine whether the su- 
persymmetric Schrodinger equations admit zero modes (they do not). For this we can 



follow Ref. 27| and notice that in a periodic potential the wave-functions are the Bloch- 



Floquet states, which for a zero-mode, if it were to exist, would translate into the condition 
4'o,±{x + L) = e^'^(j)o^±{x), with u = W{x) = (3^. Since, the real factor u is non-zero (apart 
from the state with = 0), there are no zero modes, the Witten index is zero, and therefore 
the supersymmetry is broken for our conventional s-wave superconductor. Admittedly, the 
significance of this fact is unclear (at least for topologically trivial superconductors studied 
here), but what however may be important is the fact that the existence of analytic solutions 
to the underlying Bogoliubov-de Gennes equations (16. ip should be related to the existence 
of solvable supersymmetric potentials and vice versa. Another interesting approach that 
potentially may lead to progress would be to study quasiclassical solutions to Eqs. (I6.13p . 
where the WKB method is known to work very well (it is exact in many notable cases). We 
shall however leave these questions for future work and explore below another means to treat 
the Bogoliubov-de Gennes equations to calculate the functional determinant of interest. 



C. Derivation of a Compact Expression for the Functional Determinant 

1. A Solvable Case with Non-Trivial Quantum Dynamics of the Order Parameter 

Up to this point, we have considered general properties of the Bogoliubov-de Gennes 
equations (16. ip and found just one explicit solution corresponding to the "trivial" case of 
a constant order parameter, thereby recovering classic BCS theory in this language. For 
a further progress, it is desirable to examine the properties of some other exact solutions 
in less trivial cases, but as noted above the number of known solvable cases is quite lim- 
ited. Fortunately, some additional insight comes from recent progress in the closely-related 
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problem of non-equilibrium BCS superconductivity. There is a whole class of new solutions 
that have been recently obtained that not only admit exact analytic treatment of pseudospin 
real-time dynamics, but also amazingly satisfy the mean-field self-consistency constraint (for 
some specific real-time dynamics). Even though, as we shall see below, these solutions for 
Ao('r) are not at all optimal for minimizing the imaginary-time action in equilibrium, let 
us nevertheless examine some associated exact solutions for the density matrix. We will 
present below the simplest such solution, which is the imaginary-time version of the Ansatz 
proposed by Levitov et ai— That is, let us seek the function, i?+(r), [see, Eqs. (16.71) and 
(16. 8p ] in the form 

i?+(r) = 2e/(r)-/(r), (6.14) 



where we, following Ref. 10|, identify /(r) = Aq ^(r) so that the equation for / reads 

// = /'-l- (6.15) 

This yields 

Wr) = f{r) = CO-' cos [u (r - Tq)], (6.16) 

where u and tq are arbitrary constants for the purpose of satisfying Eq. (I6.15p . However, 
we also have to satisfy the periodicity requirement for the order parameter, Ao(0) = Ao(/3), 
which leads to two possibilities: (i) If u; = 27m//3, with n G Z, the solution (I6.16P is 
"naturally periodic" with the period commensurate to /3; (ii) If 7^ 2?™//?, but tq = (3/2, it 
is an "accidentally periodic" solution. As we shall see below for the purpose of minimizing 
the imaginary-time action the latter "accidental" periodicity is much preferable, while the 
naturally periodic solution is not even allowed in the case of (I6.16p . 

In either case, the Ansatz (I6.14p immediately leads to the following solution for R+{t) = 
v{t)/u{t) 

2^ 

i?+(r) = — cos [u (r - tq)] + sin [u (r - tq)] (6-17) 

CO 

and for the function p{t) = u{t)v{t) [see, Eq. (16. 9p ] 

p(r) = 1 + ^/(20 tan [c.(r -ro)] ^,,(_„)_ 
cos [uj (r - ro)J 

These solutions (16.170 and (16.180 together with Eqs. (16. 4p . (16. 5p . and (16.70 determine the 
full "density matrix" {i.e., the solution to the original equation (16. ip as follows 



P(r) 



{e + tan(u;ro)} {e + tan [uo (r - r,)]} ; ^ ^ 

£5^; cos-Hc^ro) cos-i [uo {r-r,)]-)l + e 



2 



' "^'Ztrr'^ -^ {e-tan(c.ro)}{e-tan[a;(r-ro)]}; j l + e^'^ • ^ 
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where we introduced e = One can explicitly verify that p(r) given by Eq. fl6.19p above 

indeed satisfies Eq. (16.11) together with the initial condition /)(0) = 1. 

The solution fl6.19l) looks complicated, but for the purpose of calculating the functional 
determinant (or equivalently the "partition function," zi[Ao{t)]), we do not need its full 
form but need just its trace at r = /3. Calculating this trace, we find an interesting result 
for this particular choice of the order parameter 



OJ 



COS (cut — CUTq) 



Trp(/3) = 2cosh(e^r), (6.20) 



which as we see does not depend on Ao(t) at all (neither on frequency nor on tq) 
and is equivalent to a pseudospin not subject to any time-dependent Ao(r): /. e., 
zi[u / cos {co {t — Tq)}] = Zi[0]. This is a very curious result indeed, because it suggests 
that the functional determinant may have a much simpler form than the actual "density 
matrix" used tool to calculate it. 

To complete the analysis of this non-trivial fiuctuation, let us examine the action 
evaluated for this particular "trajectory" of Aq{t), which has the form: S'[Ao(r)] = 

^ / cos^ (^r-u^T ) ~ ^;[0]> with zi[0] given by (lOOj) . We notice that the first term di- 

verges for the trajectories with "natural periodicity" because Ao(r) changes sign, which is 
not allowed if the phase fiuctuations have been eliminated, but in any case such trajectories 
do not contribute to the partition function at all. If however tq = /3/2 and < vr (i.e., the 
order parameter is positive, Vr G [0,/3]), we find immediately that the contribution to the 
action is 



S 



CO 



^tan(w/3) + ^FG, (6.21) 



cos {lot — uj/3/2) \ g 

where the last term is the action of a non- interacting Fermi gas [c.f., Eq. f l6.20p ] and the 
first one is the energy cost to have a fiuctuating order parameter in the form (I6.16p . Since 
there is no means to compensate this energy cost by adjusting the "negative-energy" term 
associated with zi, we conclude that the chosen "trajectory" of Ao(r) is a low-probability 
event in a low-temperature superconducting state. Note that since co must be kept smaller 
than ttT, instantons of this type will completely die out in the ground state, but they may 
appear as classical excitations at higher temperatures, including even the normal state. 



2. Functional Determinant on a Class of Elliptic Functions 

Sec. IVI C II shows that while a tour-de-force derivation of the density matrix for a non- 
trivial fiuctuating order parameter is not impossible, but generally is quite complicated, 
the actual result for the functional determinant may look very simple. To understand the 
origin of the "mysterious" simplification of the complicated density matrix f l6.19p to the very 
simple- looking trace ( I6.20p , we will rely on the recent work of Yuzbashyan^^ and a related 
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work of Dzero and the author.— 

Let us consider a particular quantum trajectory of the order parameter Ao(t) and an- 
alytically continue this function from r G M to complex values z = t + it E C. One can 
formulate a sensible general framework in terms of a ^-dependent iS-matrix, S{z), but we 
will consider here only its analytically continued form on the real-time axis (which is equiv- 
alent to a Feynman-Wick rotation at T = 0). Let us now use this analytical continuation to 
relate the Bogoliubov-de Gennes equation (16.11) for the density matrix in imaginary time to 
the corresponding Schrodinger equation for the S'-matrix in real time t: 

= k{it)Si{t) = ^^^^^^^ j m, with 5,(0) = i. (6.22) 

where Ao(zt) and hi{it) are symbolic notations for the dynamic order parameter and the 
Hamiltonian properly analytically-continued to real times, correspondingly. The Hamilto- 
nian, hi{it) = Re Ao(it)r^ + Im Ao(«t)'r^ + ^^r^ = (b, ■ f) /2, belongs to the two-dimensional 
representation of the Lie algebra, su(2), while the unitary S'-matrix belongs to the two- 
dimensional representation of the SU{2) group. Note that there is no need to specify the 
dimensionality of the matrix representation of operators in Eq. fl6.22p . which can be viewed as 
an equation of motion in the abstract group, i.e. hi{it) G su(2) ~ so (3) and S{t) G SU (2), Vt. 
We can also write an associated Schrodinger equation for spinor wave- function, ^ = (J^ 

Just like in Sec. IVI At we can argue that if we know a particular solution to Eq. fl6.23p that 
satisfies an arbitrary initial condition, we can always construct another linearly independent 
solution with the help of a time-reversal operation that now reads (— zf^) \E'*(t) (zf^) and 
hence a full S- matrix can be constructed using one particular solution to (16.231) . Let us note 
here that a key motivation for studying the analytically-continued form of the Bogoliubov-de 
Gennes Eqs. (16. ip and (16.20 is that if we know the solution to Eq. (16.220 or just a particular 
solution to Eq. (I6.23p . which has the familiar form of a Schrodinger equation, we should 
be able to analytically continue the result back to imaginary time [so that in some sense 
S{—iT) — > p(r)] and therefore calculate the functional determinant. 

Now, let us narrow down a class of functions considered from that of arbitrary order 
parameters (analytically continued from the values r G [0, f3] to complex arguments), Ao(z), 
to those that are periodic along the it-axis, i.e., Ao(-2) = Ao(z + i/3t), where /3( G M is the 
corresponding period. Let us also assume that the original fluctuation is a meromorphic and 
periodic function along the r-axis as well, i.e., Ao{z) = Ao(-2 + (3r), which is either due to 
a natural periodicity with the period commensurate to the inverse temperature /3t- = (3/n 
or some other period unrelated to the fact that the relation Ao(0) = Ao(/3) must hold (see. 
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e.g., previous Sec. IVI C H where depending on the parameters both cases can be reahzed). 
This narrower class of functions represents elhptic functions,— with primitive periods and 
iPti which are arbitrary constants at this point. 

If ^o{z) is an elhptic function with the primitive periods (PrjiPt) as defined above, 
then the Schrodinger equation (16.221) is that describing a spin-1/2 under a periodic-in- 
time perturbation and so let us look for a solution of (I6.22p in a Bloch-Floquet-type form, 
S{t) = Sp{t)e'^^ + Sl{t)e-'^\ where Sp{t) = Sp{t + A) is a periodic 2 x 2-matrix and £ is a 
constant. On the other hand, we could have used the same Floquet argument for the original 
Bogoliubov-de Gennes equation to argue that the "density matrix" may be written in a simi- 
lar form, /3(r) = pp +(r)e^'^+Pp _(r)e~^^, where now Pp,±(t) = Pp,±{t+(3r) is a "periodic part 
of the density matrix" [c.f., Eq. (I6.19p ]. and S is the same as before. These arguments sug- 
gest themselves to be generalized in the form of a solution <S{z) = <Sp^^{z)e^^ + Sp-{z)e~^'' , 
where z = r + it and Sp^±{z) is an "elliptic matrix function" of a complex argument, z G C, 
such that <Sp^±{z) = <Sp^±{z + nf3T + iml3t), Vn, m G Z. 

If we now assume that the periodicity of Aq{z) along the r-axis is commensurate with 
the "natural" period, /3, then we find immediately that Sp^±{z) = 1/2 because of the initial 
condition 5(0) = /)(0) = 1 and the imposed periodicity and therefore we conclude that 



where the current assumptions are that Ao(z) is an elliptic function with the periods (/?,-, i(3t) 
such that f3//3r G Z and f3t is arbitrary and the solution of the Schrodinger equation (I6.22p 
has a Bloch-Floquet form. 

To get a more useful expression for the "partition function," z, let us now focus on a 
dependence Ao(it) slow enough so that there are no level crossings taking place. Consider 
now a pseudospin, described by the spinor \E'i(t) of Eq. (I6.23p . evolving from the initial 
state that is an eigenstate of the Hamiltonian at t = 0, e.g., we can take it to represent to 
a pseudospin moment opposite to the "initial magnetic field," hi{0) = (Ao(0), 0, ,^), which 
lies in the XZ-plane. The adiabaticity assumption immediately tells us that the quantum- 
mechanical phase "collected" after the completion of a single cycle, t : — )■ /3t of the 
magnetic field, b(t), is given by the expression in the exponential below 



where 7Berry is the Berry phase, which is determined by the following flux through the area, 
swept by b(t) over the cycle, t : — )■ 



z[Ao(r)] = 2cosh(^/3), 



(6.24) 



^i(/3^) = ^i(o)e~*(^«^"^+^^^°). 



(6.25) 




(6.26) 
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and 7dyn is the dynamical phase given by 



7dyn 



e + \Aom'dt, 



(6.27) 



where the integrand can be easily recognized as an "instantaneous eigenenergy" of the cor- 
responding spin Hamiltonian (which in turn represents the energy of an excitation in a 
superconductor subject to such a fluctuation). Note that we could have taken the other 
initial condition corresponding to a pseudospin pointing along the "initial magnetic field," 
which would have evolved into a state with the dynamical phase, which is a complex con- 
jugate to (z7dyn) above. We can recall now that since we are interested in the S'-matrix 
modulo its periodic part, we can construct the remainder out of the two phase factors, 
which therefore gives exactly the desired {S/3t) that appears in Eq. fl6.24p . 

If we now make a further simplifying assumption and consider a fluctuation, Ao(t — tq), 
which is described by an even function of its argument (see, previous Sec. IVI C ip such that 
the analytically-continued Ao(zt) is also real- valued, we immediately find that the effective 
"magnetic field" simplifies to hi{t) = (Ao(it), 0, ^). Therefore, the area swept by any such 
dependence in the parameter space is zero and the Berry phase fl6.26p vanishes identically 
as well. Note that this conclusion would also hold if we assume that the order parameter 
is an odd function of (r — tq), such that it leads to a purely- imaginary Ao{it) (let us recall 
that phase fluctuations in imaginary time have been eliminated). Under these assumptions, 
we can identify the factor arising from the "non-periodic" part of the S'-matrix/ "density 
matrix" with the dynamical phase, {Sf3t) = 7dyn to obtain the following result for the 
functional determinant: 



zi[Ao{t)] = 2 cosh 



(6.28) 



where we have restored the index / that parameterizes the sites of the original Richardson 
model. 



3. The Adiabaticity Requirement 

Eq. fl6.28p has appeared after a chain of rather general arguments, which however included 
a number of additional assumptions. Let us reiterate these assumptions: We have assumed 
that Ao(t) can be analytically-continued from r G [0, /3] G M to the complex plane, C, and 
that the resulting function Ao(-2) is an elliptic function with two primitive periods {I3r,i(3t) 
such that f3/f3r G Z. We also have further assumed that the there exists tq G M such that 
the order parameter is either an odd or an even function of the argument A(r — tq), which 
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ensures that the Berry phase vanishes. Finally, we have assumed that the time- dependence 
of the analytically-continued order parameter, Ao(zt) is "slow enough," such that no level 
crossings take place. 

The last assumption is the most restrictive and one may wonder about the accuracy 
and domain of applicability of the conjecture (I6.28P and in particular about the meaning of 
"slow enough" in the adiabaticity assumption. This question has been addressed by Gan- 



gopadhyay, Dzero, and the author in Ref. |25| in the context of two-level-system dynamics 



in superconducting qubits. Mathematically, Ref. |25| presented an extended class of exact 
solutions associated with elliptic functions describing the driving field (which represent a 
generalization of the anomalous solitons discussed in the amazing paper of Yuzbashyan in 
Ref. Here we reiterate only key facts relevant to our paper: 

The following functional dependencies of Ao(^^) admit exact explicit solutions of the 
associated Eqs. flOSD and fl6:22|) : 

A / N ^ ^ 1 — ?7+sn^ (wt, k) , , 

A„(e..i) = n. + n. ^^:;;^^,;^^;^j , (6.29) 

where e = (ei, 62, 63) is a shorthand to describe three parameters that appear in the following 
equation [c./., Eq. fl6.15p ]. P = {f — ei)(/ — e2)(/ — 63), and are subject to the constraint 
ei + 62 + 63 = 0. A solution to the equation for / above can be expressed in terms of 
the Weierstrass elliptic function, which is related to the order parameter f l6.29p . The other 
constants in Eq. fl6.29p : Qa, ^, and k, are not free but are determined uniquely by 
[E.g., K = (e2 — e3)/(ei — 63)]. Finally, the function sn (■,«;) in Eq. f l6.29p is the doubly- 
periodic Jacobian elliptic function. 

The derivation of the associated solution involves Bloch representation of Anderson pseu- 
dospins as follows M{t) = l-^^{t)T-^{t), where of course M'^{t) = 1/4 so that M{t) G S"^, 
which is the standard Bloch sphere. The equations of motions for the Bloch vectors follow 
from the Schrodinger equation in a standard way and yield the familiar Bloch equations: 

M{t) = -z^t(t) [f , hi{it)] ^{t) = hi{t) X M{t) (6.30) 

where the "magnetic field" in the cross product is exactly as in Sec. IVI C 21 hi(t) = 
(Ao(it), 0, (^i) with Ao(it) given by Eq. (16.291) . Therefore, the equations of motion (I6.23P 
in the 5'f/(2)-group have been reduced to equations of motion (I6.30p on a sphere, 5*^. Let us 
recall that S"^ = SU{2)/U{1), therefore Eq. (I6.30p has less direct information than the orig- 
inal Schrodinger equation. It turns out that the "missing part" is exactly the sought-after 
overall time-dependent U{1) phase of the wave-function, which we expect to reduce to the 
sum of the Berry phase and a dynamic phase discussed in the previous Sec. IVI C 21 [in the 
case of dependence (16.290 . the Berry phase is zero]. 



Using the Ansatz proposed by Yuzbashyan in Ref. [12], one can find^ the solutions. 
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M(t), to Eqs. f l6.30p . expressed in terms of elliptic functions with the same periodicities that 
the elliptic function fl6.29p in accordance with the suggested generalization of the Floquet 
argument to eUiptic functions, as discussed in the previous Section IVI C 2[ Using these 
exact solutions one can construct the full S-matrix describing the motion in SU{2). It can 
be done by parameterizing the components of the spinor in Eq. f l6.23p as follows ip^/_i{t) = 
|ctt/4.(i)| e*'''^*''^^^*^*^. One can see that while the amplitudes and relative phase are directly 
related to the "instantaneous" direction of the Bloch vector \a-^/i{t)\ = ^l±2M'{t) and 
6{t) = arctan [M^(t)/M^(t)], the common phase 7(t) depends on the trajectory in a non-local 
way and to determine it, one has to go back to the Schrodinger equation (16.230 . This indeed 
can be done and the phase, 7, can be found (this part has to be done numerically for generic 
parameters). The main conclusion of this analysis is that if Qa is small, this exact phase 
is essentially indistinguishable from the dynamic phase described by Eq. (I6.27P [however, 
for any non-zero fi^, Eq. (I6.28P is not exact]. As increases up to a critical value, l^i"'', 
level crossings start taking place and the absolute value of the quantal phase is suppressed 
compared to the adiabatic result. In all cases considered, the adiabatic quantal phase is 
either equal or larger than the exact phase, and therefore it can be viewed as an estimate 
from above. Ref. 25] also indicates that for all ila < ^'a^\ the adiabaticity condition is 
satisfied and in this case the compact (16.280 expression for the functional determinant can 
be used. 



4- Explanation of the "Fermi Gas" Result Obtained in the Solvable Case of Sec. VVI C li 



This discussion was initially motivated by the "paradox" found in the fully-solvable case 
described in Sec. IVI C II We remind that the contribution to the action of a pseudospin 
moving in the presence of a non-trivial fluctuating Ao(t) described by Eq. (16.160 turned 
out to be completely independent of the parameters of this fluctuation and was found to be 
identical to the corresponding contribution expected in a Fermi gas, i.e., in the absence of 
any order parameter whatsoever, A(r) = 0. Another part of the "paradox" was that the full 
solution for the "density matrix," p(t), was very cumbersome (I6.19P and the simplification 
occurred at the final stage of calculating its trace, which led us to Eq. (I6.20p for Trp;(/3) = 
2 cosh {C,it). 

These paradoxes can be now resolved with the help of Eqs. (16.281) and (16.290 . One can 
consider various limits of the function (I6.29p . in particular, that of — 0, which leads to the 
following expression of the dynamic order parameter considered previously by Yuzbashyan^^ 
and Levitov et al.-M 

= w dn Lot, 



Ao(e,it) 



(6.31) 



On the other hand, the limit k — )■ 1 leads to sn {u, 1) = tanhw and Eq. (I6.29P reproduces 



the anomalous soliton of Ref. 



12 



If we now take both limits, 



I.e.. 



K — )■ 1 and ^la — )• 0, we 
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find 

Ao{e,tt) = ^ (6.32) 

An analytical continuation of this function to imaginary time yields Ao(e,r) = 

oj/ cos(a;t), which is exactly the soliton studied in Sec. IVI C II Note that this soliton is not 
an elliptic function, but a circular function, because it has only one incommensurate period 
in the r- "direction." However, it is a limiting case of a proper elliptic function, with its 
period (it taken to infinity. Therefore, per the arguments of Sec. IVI and using Eq. fl6.28p . 
we are to write the "partition function" as follows: 



zi[uj/ cos {ojt)] = 2 cosh 




2 cosh (/3^;)- (6-33) 



Now, it is easy to see the origin of the paradoxical result fl6.20p : Since cosh^^(a;t) decays 
exponentially with increasing t, the soliton term in Eq. fl6.33p vanishes in the {Pt — > oo)-limit 
and does not contribute to the integral. We therefore recover the correct result fl6.20p ! 

Now that the origin of the result (16.200 in the exactly-solvable case is understood, we can 
use the exact solution to get another insight into the range of applicability of Eq. (I6.28p . 
The list of assumptions for the validity of Eq. (I6.28P includes that /3 be a natural rather than 
an accidental period of Ao(t) (enforcing this periodicity via a periodic repetition of Aq{t) 
from r G [0, /3] — i- M would not necessarily work, because the resulting periodically-continued 
function may not have the desired analytical properties and any arguments based on them 
would become unreliable). Hence, we do not expect Eq. (I6.28P to work in the accidentally- 
periodic cases, but the exactly-solvable example (I6.16P with tq = [3/2 and Wu G M shows 
that at least in this particular case Eq. (16.281) does work correctly. It remains unclear at this 
stage, whether this result is an artefact of the particular dependence (I6.16P or it is rather an 
indication that Eq. (I6.28P applies to a wider class of elliptic functions and their limits with 
accidental /3-periodicity. 

D. Contribution of Elliptic Trajectories to the Partition Function 

Now let us summarize our findings (conjectures) and present the following expression 
for the contribution to the partition function of those specific elliptic trajectories, Ao(r) — >■ 



28 



Aq{z) G Ell, for which our expression for the functional determinant applies: 



-/3 E 

oZeu = e 



V 



Ao(2)eEii 



S[A{t)]>^ I \Ao[r) 





"A|(l)" 
_ 27rg 

2 5: In 



g-5[A(r)] ^j^g^g 



(6.34) 



cosh 



2A 



+ |Ao(^t)|'c?t + ^7Bcrry [Ao(^t)] 



where just as before Ao(it) corresponds to an analytically-continued order parameter, which 
leads to an elliptic function with the primitive periods {(3r,iPt) or a limit of such an elliptic 
function. In Eq. fl6.34p . we have also included the Berry phase contribution, which in general 
should be present, but whenever Ao(it) is either purely real or purely imaginary, the Berry 
phase vanishes identically. 

We have already verified that the action in Eq. fl6.34p reproduces the exact results in 
certain exactly-solvable cases. Note here that the classic BCS result (14. 9 p is certainly re- 
produced exactly as well, because a constant function represents a trivial elliptic function, 
Ao(r) = Ao(it) = A = const and so we can use Eq. (I6.34p . and after an integration recover 
the correct "partition function" of a spin in a constant mag netic field, = 2cosh(E//3), 
which corresponds to the BCS mean-field. One can also argue based on (I6.34p that the 
classical mean-field is indeed a true minimum on the space of these elliptic functions. Ell. 
Let us consider an order parameter Ao(r) = A + 5A(r), where SA{t) is in some sense small. 
Let us expand the action in Eq. (I6.34p assuming that 6A{t) does not induce a Berry phase. 
Hence a correction to the relevant part of the action is (we consider the low-temperature 
limit, (3 — 7- oo): 



SSo 



lac 



fit 



-At 



(6.35) 



it + Ao 



Note now that the integral above can be equivalently written as an average of 5A{it) over an 

/3t nPt 

infinite number of periods, J 5A{it)dt = lim ^ J 6A{it)dt. Therefore, the contour of 
T+it is that going from to ioo. The difference between this integral and 



integration over z 
that going along the r-axis is 
i.e., it is given by the sum of al 



;-2vrz)E..res5A(zO = 0, 



j5A{it)dt- J 5A{T)dT 



residues of 6A{z) enclosed in the first quadrant. It is equal 
to zero, per one of the elementary properties of elliptic functions,^- which states that the 
sum of all the residues of an elliptic function inside a period-parallelogram always vanishes. 
Therefore, one can write the first variation of the both parts of the action in terms of the 

(l/^)-AoEfe + Ar"'^'^ 



same function, 5A{t): 5S = 2Aq 
is satisfied for the BCS mean-field. 



Z6£ 
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One may wonder, if one can use the analytical properties of elliptic functions to bring 
the integral that appears in Eq. fl6.34p to the r-axis in a similar way, i.e., to the form 



/3 I 

— )■ J y + \ Ao{t)\ 7 We know however that this substitution can not generally be cor- 



rect, because some non-trivial solutions that we have analyzed manifestly contradict this 
assumption. However, we have proven above that for all relevant fluctuations in the 
immediate "vicinity" of the classic BCS mean-field the substitution above would work. 
One can explicitly verify that the interesting property of this (generally incorrect) sub- 

stitution is that the variational analysis of the functional, S* [/o(t)] = / '^t||/o(t)| — 





lec 



cosh| i/^/ef + |/o(r)rrfr 



[i.e., the constraint |^ = 0) indeed immediately 



selects the classical mean- field /o(t) = Abcs = const as the only saddle point. Hence, one 
can use the expression for S above to determine the contributions to the partition function 
due to Gaussian quantum fluctuations in the vicinity of the BCS mean-field (for simplicity, 
we consider the low-temperature limit only). The result is not unexpected and is quite 
"boring," taking the following form for the usual BCS superconductor {i.e., the parameter 
space, £, is momentum space): 



^nearBCS 





'5A(r)" 




ng 



/3 

■ fSA^Mdr 

, (6.36) 



where J^bcs is the energy of the classical mean- field BCS state, T = /3~^, and V is the actual 
physical volume and hence the contribution of these mesoscopic fluctuations to observables 
in a bulk system is negligible. 

It is alluring to attempt a variational analysis of the action in Eq. (16.341) to see if 
there could exist other saddle points apart from the classical mean-field. However, the 
variational analysis would be problematic, because the action in Eq. f l6.34p contains "ap- 
ples and oranges," that is two functionals of different functions, Ao(t) and Ao{it), which 
are related in a non-trivial way via an analytical continuation. However, to make the 
case that non-linear soliton contributions are important (we should distinguish here be- 
tween instantons, which are trajectories that connect classical minima,— and fictitious at 
this stage new minima, which we dub solitons), one does not necessarily need to find 
true quantum minima, finding any quantum trajectory that corresponds to the energy 
smaller than mean-field would suffice. To clarify the content of this (open) problem, let 

US introduce two parameters = J Al{T)dT and A2 such that ( f + Ag ) = 

^ ^ 

^ J \J + \Ao(it)fdt + ^'jBerry [Aq (it)] j . In thcsc notations, the action in Eq. fl6.34p 
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takes the form [c./., Eq. (14. 9 p 



5[A„A2] = ^-2Vln 

Since the first term is always positive and the second one is always negative, we are to look for 
ways to minimize Ai and maximize A2. In the classical BCS mean-field Ai = A2 and there is 
no room for any additional variation, but in the functional f l6.34p such additional variations 
are in principle allowed. One can check that the analytical continuation of solitons of the 
non-equilibrium BCS problem with natural periodicity do not produce a "good" solution at 
least at T = 0. If on the other hand, we take Eq. (16.341) as a given functional and consider 
various trial functions without attempting to prove that they actually satisfy the formal 
domain of validity of the Ansatz, we immediately find a variety of dependencies that do 
"better than classical mean-field" in term of energetics. However, these "results" should be 
taken with a grain of salt, because there is no way to determine the actual range of validity 
of ( 16.34P beyond those dependencies associated with known integrable spin dynamics and 
the most natural explanation for any accidental solution obtained within a trial-and-error 
analysis of (I6.34p is that it is probably beyond the applicability of the method. On the other 
hand, there appears to exist no proof that such solitons are impossible. Looking at the rich 
structure of the functional determinant, it appears conceivable that there exist trajectories 
in the huge functional space spanned by, A(r), that do not just collapse into the mesoscopic 
term (16.361) . but instead provide more noticeable contributions to the action. A numerical 
analysis of some non-linear solutions, guided by the analytical result (16.341) . will be published 
elsewhere. 

VII. SUMMARY 

This paper presents an analysis of non-perturbative fluctuation phenomena in the pairing 
model. The key step of this analysis is a decomposition of the partition function of the 
Richardson model into spin and pseudospin terms. It is shown that such factorization is 
possible for a generalized Richardson model that includes both BCS and spin interactions. 
Even though we have not presented here a theory to describe both types of non-trivial 
interactions on an equal footing, the development of such an extension is straightforwarci^ 
and would lead to a two-order-parameter theory expressed in terms of two "global" Hubbard- 
Stratonovich fields.— The analysis of phase fluctuations presented here indicates that these 
interactions will be competing and that such competition can be enforced via a commutation 
relation between the density and spin density and the overall phase. However, the present 
paper has focused on the analysis of a simpler canonical Richardson model that has no 
magnetic interactions. Even though the spin sector of this Richardson model is trivial, the 
existence of this (single-particle) sector is important for the possible existence of any non- 



2 cosh 



V^fTAf/? 



(6.37) 
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trivial fluctuations of the amplitude of the order parameter in the low-temperature phase. 

The main technical part of the paper involves a calculation of functional determinants 
that appear in the non-linear effective action expressed in terms of the Hubbard- Stratonovich 
field. We have shown that the Anderson pseudospin language and in particular its coherent- 
state path-integral representation lead to practically useful and physically intuitive insights 
into the structure of the functional determinants for non-trivial quantum trajectories. There- 
fore, this approach may be much preferable to the conventional Grassmann path integral 
method. We have shown that a functional determinant is given by the trace of a density 
matrix that satisfies the Bogoliubov-de Gennes equations in imaginary time. This leads to 
a differential equation of the Riccati type, which is directly related to the supersymmet- 
ric Schrodinger equation with superpotentials determined by the imaginary-time dynamics 
of the order parameter, A(r). Let us note here that a particularly promising direction for 
further research could be to use the WKB-method to treat the relevant differential equations. 

In Sees. IVI C 21 and IVI Dl we proposed an explicit, compact expression for the functional 
determinant for a certain large class of elliptic functions and the arguments that led us 
to the conjecture (16.341) involved an analytical continuation of the Bogoliubov-de Gennes 
equations in imaginary time to the real-time axis (or more generally to the complex plane, 
z = T + it), such that the problem could be mapped onto that of a two- level system in a 
time- dependent magnetic field determined by quantum dynamics. This is a known, very 
complicated problem, but we have taken advantage of some recent exact results and our 
recent work on an extension of these results to analyze a family of exact solutions that are 
associated with elliptic functions. These results have led us to Eq. (16.341) . which provides 
a useful intuition for the effective action of the model and suggests that the functional 
determinant, that is often treated as a thing- in- itself, can actually be calculated and is 
related in a very straightforward way to the dynamical and Berry phases of a pseudospin 
moving in a "magnetic field," determined by the quantum dynamics of a fluctuation. Let 
us reiterate however that a formal justification of our solution applies only to adiabatic 
dependencies on the specific class of elliptic functions with the periods along the imaginary 
and real-time axes. 

An important open question is whether the considerations presented in this paper can 
be generalized to other types of functions, A(r), which are not associated with any elliptic 
functions that lead to integrable pseudospin dynamics. A particularly promising avenue 
here could be to use the reverse-engineering approach for constructing exact solutions as 
described in Refs. and 30 1, which effectively implies a change-of- variables from the 
Hubbard-Stratonovich field, A(r), to the generators, $(t), which govern the dynamics of the 
S'-matrix, Sit) = exp [— |$(t) ■ <t] , satisfying the proper Bogoliubov-de Gennes equations. 
It would also be interesting to see whether a chaotic rather than integrable dynamics^ can 
be realized under any circumstances in this model. Quite generally such dynamics, if at 
all possible, are not expected to lead to energetically favorable contributions to the action, 
because the "trivial" term that gives an energy penalty to any non-zero order parameter 



32 



configuration corresponds to the average of |Ap, wliile tlie second non-trivial term tfiat 
favors superconductivity contains contributions from different sites, and if tlie dynamics 
exliibit a "cliaotic beliavior" in tfie parameter space, £, tlie signs in the second term would 
fluctuate strongly from site to site and are expected to average out to zero instead of lowering 
the corresponding energy. This argument supports the approach to use regular elliptic 
trajectories that describe a synchronized collective behavior of the pseudospins. Another 
open question relates to the role of the Berry phase in the functional determinant f l6.34l) . All 
exact solutions we have analyzed [that are sensible to describe thermodynamics, where the 
constraint A(0) = A(/3) must be imposed] have a trivial (zero) pseudospin Berry phase. This 
however represents a limitation in our ability to solve Eqs. f l6.ip . rather than an indication 
that Berry phase terms are unimportant. 

Finally, we reiterate the main question posed in this paper and the arguments of the last 
Sec. IVI Dl which suggest that non-perturbative soliton trajectories that co-exist with the 
classical mean-field are not impossible and in fact the rich general structure of the functional 
determinant suggests that the construction of such quantum fluctuations may be possible 
at least in some modification of the model (which may involve interactions for real spins). 
Generally, the right question to ask would be whether there exists any fermion model that 
exhibits breaking of continuous symmetry and such that its low-temperature phase allows 
non-perturbative soliton solutions for a component of the Hubbard-Stratonovich field that 
is normally considered "massive?" In other words, can the non-linear effective action for 
the Hubbard-Stratonovich field develop any other minima apart from the classical mean- 
field? A proof that no such solutions exist would confirm the fundamentals of classical 
spontaneous symmetry breaking and mathematically would imply that there is no need 
to study complicated non-linear actions at T = [such as the non-linear effective action 
in Eq. (12. 6p ] and that in an infinite system they should crossover to a functional delta- 
function of the type, e-^t^Wl oc 5 [|A(r)| - Amf] , c./., Eq. (I06|) . On the other hand, 
even a single example of an order-parameter trajectory that is energetically beneficial to the 
classical mean-field would seriously question this fundamental conjecture. We know that 
any such trajectory, if at all possible, can not be anywhere near classical mean- field (in the 
functional space of allowed fluctuations), but the possibility of a non-perturbative solution 
not adiabatically-connected to the mean- field has certainly not been ruled out. 
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